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Small Bites: Star formation recipes in extreme dwarfs 
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ABSTRACT 

We study the relationship between the gas column density (Shi) and the star formation 
rate surface density (Esfr) f° r a sample of extremely small (Mb ~ —13, AV50 ~ 30km s _1 ) 
dwarf irregular galaxies. We find a clear stochasticity in the relation between the gas column 
density and star formation. All gas with Ehi > 10M Q pc -2 has some ongoing star formation, 
but the fraction of gas with ongoing star formation decreases as the gas column density de- 
creases, and falls to about 50% at Shi ~ 3M© pc~ 2 . Further, even for the most dense gas, the 
star formation efficiency is at least a factor of ~ 2 smaller than typical of star forming regions 
in spirals. We also find that the ratio of Ha emission to FUV emission increases with increas- 
ing gas column density. This is unlikely to be due to increasing dust extinction because the 
required dust to gas ratios are too high. We suggest instead that this correlation arises because 
massive (i.e. Ha producing) stars are formed preferentially in regions with high gas density. 
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O 1 INTRODUCTION 

1 ■ 

Models of galaxy formation and evolution generally use semi- 
empirical "recipes" to follow the process of star formation (e.g. 
Springel et al. 2005; Governato et al.l l2010). Typically, star forma- 
tion is assumed to set in only above a "threshold" gas (column) 
density E gas and beyond that to be proportional to a power of E gas . 
This is supported by observations of nearby star forming galaxies 
(e.g. ISchmidlll 19591 ; lK^nnicutllll998t) . However, most of these ob- 
servations are of large spiral galaxies, whereas from the hierarchi- 
cal galaxy formation model one would expect that the first formed 
systems were much smaller than the typical z ~ spiral. Here we 
study the relation between gas and star formation in nearby, ex- 
tremely faint ( Mb ~ —13, AV50 ~ 30 km s" 1 ) gas rich dwarfs. 

The dwarf galaxies in our sample are dynamically and struc- 
turally very different from the large spiral galaxies for which the 
widely used star formation recipes have been derived. Firstly, in 
our sample galaxies the rotati on velocity is not much larger than 
the velocity disp ersion (e.g. Begum, Chengalur, & Hoppj 120031 : 
lBegumetal.ll2008h . Further, the gas does not settle into a thin 
disc; the mean observed a xial ratio of the gas discs is ~ 0.6 
( Rovchowd hurv et al.l|201o| ). Both this as well as the expectation 
that negative feedback from su pernovae would play a mor e impor- 
tant role in small galaxies (e.g. lMac Low & Ferrarall 19991) make it 
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likely that the relationship between the gas density and star forma- 
tion in dwarf galaxies is different from that in spirals. 



Observationally, there is another major difference between 
studies of star formation recipes in dwarf galaxies and spi- 
rals. M olecular gas is almost never detecte d in dwarf galax- 
ies (e.g. iTavlor. K obulnickv. & Skillmanl l 19981) , which means that 
the gas column density has to be estimated from the HI col- 
umn density Ehi alone. On the other hand, in large spirals, 
the star formation appears to be governed by the molecular gas 
density and to be much les s (if at all) r elated to the atomic 
gasfe.g. IWong & Blitzl |2002| ; llxrov et al.1 120081) . However, in 
iRovchowdhurv et alj {2009) (henceforth R09) we showed that for 
dwarf galaxies, in regions of active star formation, the star forma- 
tion rate Esfr is correlated to the HI column density, albeit with 
significant scatter. R09 also found that there was no sharp "thresh- 
old" for star formation, with star formation proceeding at all gas 
colum n densities, down to the sensitivity limit of the data. Simi- 
larly, iBigiel et alj J201oh find that in the HI dominated outskirts of 
spiral galaxies, the SFR and HI are correlated, albeit with a scatter. 



In this paper we extend our previous work in two important 
directions. Firstly we try to quantify the stochastic nature of the 
relationship between Esfr and Ehi- Secondly we also study the 
relationship between the gas column density and the formation of 
stars of different masses. 
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Table 1. The sample 



Galaxy 


M B 


Dist 


Group 


a 


b/a 




(mag) 


(Mpc) 




(') 




UGC 685 


-14.31 


4.5 


Field 


1.4 


0.71 


KK 14 


-12.13 


7.2 


N672 


1.6+ 


0.37 


UGC 3755 


-14.90 


6.96 


Field 


1.7 


0.59 


KK 65 


-14.29 


7.62 


Field 


0.9+ 


0.56 


UGC 4459 


-13.37 


3.56 


M81 


1.6 


0.87 


UGC 6456 


-14.03 


4.3 


M81 


1.5 


0.53 


KK 144 


-12.59 


6.3 


CVnl 


1.5+ 


0.33 


DDO 125 


-14.16 


2.5 


CVnl 


4.3 


0.56 


UGC 7605 


-13.53 


4.43 


CVnl 


1.1 


0.73 


DDO 181 


-13.03 


3.1 


CVnl 


2.3 


0.57 


KKH 98 


-10.78 


2.5 


Field 


1.1+ 


0.55 



+ : diameters correspond to the Holmberg system ( 26.5 mag arcsec 2 ) 
2 SAMPLE AND DATA ANALYSIS 

Our sample consi sts of 23 galaxies d rawn from the GMRI0 FIGGS 
HI 21cm survey teegum et al]|2008h with UV data from GALEX^ 
See R09 for details. For 1 1 of these 23 galaxies Ha data from the 
6m BTA telescope in Russia is available. The full sample has me- 
dian HI mass Mhi ~ 28 x 1O 6 M0, median blue magnitude Mb ~ 
— 13.2, and median velocity width AV50 ~ 32 km s" 1 . The cor- 
responding values for the Ha subsample is Mhi ~ 34 x 10 6 Mq, 
Mb ~ —13.5, AV50 ~ 33 km s" 1 . The galaxies with Ha ob- 
servations are listed in Table [T] the columns in the table are: Col- 
umn(l) the galaxy name, Column (2) the absolute blue magnitude 
(corrected for galactic extinction, the internal extinction correc- 
tion has been assumed to be negligible), Column(3) the distance 
in Mpc, Column(4) the grou p membership of th e galaxy. All of 
this data has been taken from Begum et alj d2008f) . Column(5) the 
de Vaucouleurs (25 mag/arcsec 2 ) diameter of the optical disc. For 
dwarf low surface brightness galaxies from the KK lists (KK14, 
KK65, KK144, KKH98), the diameters correspond to the Holm- 
berg system ( 26.5 mag arcsec -2 ). Column(6) the optical axis ra- 
tio. Data for columns (5) and (6) have been taken from taken from 
iKarachentsev et alj {2OO4T) . 

Background corrected GALEX FUV band (1350-1750 A) 
images were converted into luminosity units using the cali- 
bration information provided at the GALEX site. Correction 
for galactic extinction was done using extinction values of 
Schlegel, Finkbeiner & Davis ( 1 19981) and using formulae from 
Cardelli. Clayton & MafhisUl989l) to extrapolate to the FUV band. 
No correction for internal extinction was made, since our sample 
galaxies are expected to be extremely dust poor. The luminosity 
values thus obtained were conver t ed to s tar formation rates using 
the calibration given in Kennicutt ( 19983) : 

SFR(M Q year" 1 ) = 1.4 x 10" 28 L„ (ergs s" 1 Hz" 1 ) (1) 

1 We thank the GMRT staff for having made possible the observations used 
in this paper. The GMRT is run by the National Centre for Radio Astro- 
physics of the Tata Institute of Fundamental Research. 

2 Some of the data presented in this report were obtained from the Multi- 
mission Archive at the Space Telescope Science Institute (MAST). STScI 
is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. Support for MAST for non-HST 
data is provided by the NASA Office of Space Science via grant NAG5- 
7584 and by other grants and contracts. 



In deriving this calibration it is assumed that the stellar distribution 
has solar metallicity and a Salpeter IMF, and that the galaxy has had 
continuous star formation over time scales of 10 8 years or longer. 
The implications of these assumptions are discussed in Section|3] 

Details of the Ha data reduc tion can be obtained from 
IKarachentsev & Kaisirj 120071) and lKaisin & Karachentse^ (20081) . 
The images were corrected for dust extinction due to our own 
Galaxy in a similar way as was done for the FUV maps. The Ha lu- 
minosit y was converted to star formation rates using the calibration 
given in lKennicutil ll998al) : 

SFR(M year" 1 ) = 7.9 x 10~ 42 L Hq (ergs s" 1 ) (2) 

The assumptions used to derive this calibration are the same as that 
used in deriving the FUV flux -SFR calibration 

For data from all the three wavelengths, relevant parameters 
(Ehi and Esfr) were calculated over several scales, viz. a) an 
average over the entire star forming disc of the respective galaxy 
(i.e. "global" values). The "star forming disc" is defined as that 
within the radius at which the star formation rate is 1.85 xl0~ 4 
M0 yr" 1 kpc" 2 (as measured from the FUV flux, with the GALEX 
images smoothed to 400 pc linear resolutions). This approximately 
corresponds to the B band Holmberg diameter for those sample 
galaxies for which the Holmberg diameter has been measured, b) 
"pixel" values. We use "pixels" that Nyquist sample squares 400 pc 
or 150 pc in size. For the HI images 400 pc resolution images are 
available for all the galaxies in our sample. Similarly for the FUV 
data, 150 pc resolution images are available for all galaxies. 

Figure Q] shows Ha greyscale images overlayed with FUV 
and HI contours, for a representative galaxy in our sample. 



3 RESULTS AND DISCUSSION 

Figure [2jA] shows the relationship between the disc-averaged 
Espjj, (and corresponding Esp^) and Ehi for the galaxies in our 
sample with Ha data. Note that the galaxies are forming stars 
even though their typical gas density is at or below the "thresh- 
old density". Panel [B] shows how the global Esfr estimates 
obtained using the two different tracers relate. Although the SFR 
tracers do correlate, there is a considerable scatter about the 1:1 
line. Note th at the data agrees b etter with the original calibration 
suggested by Kennicu tt| d 1998a ) than with the re-calibration sug- 
gested bv lLee et elj 120091) . though it should be noted that the latter 
sample is much larger than ours. In terms of total SFR the values 
range from 2.79 x 10" 4 M yr" 1 to 1.05 x 10" 2 M Q yr" 1 
with FUV as tracer, and from 2.33 x 10" 4 M yr" 1 to 
1.21 x 10~ 2 M yr" 1 with H a as tracer. Finally, following 
iHunter. Elmegreen & Ludkal d2010h we show in Panel [C] the ratio 
I^fr'/Esfr as a function of Eg^. There is a clear correlation 
and the best fit line has a slope of — 0.63±0.09, (compared to -0.59 
obtained by Hunter, Elmegreen & Ludka (2010)). The SFR calibra- 
tion w e used assumes solar metallici t y, how ever, as discussed in de- 
tail by Hunt er. Elmegreen & Ludka (2010), the fact that the dwarf 
galaxies have lower than solar metallicity has only a marginal ef- 
fect on the Eg^jJVEgjJjj ratio, since both calibrations are similarly 
affected. 

In what follows we take a look at the relationship between 
gas and star formation on small scales, by making "pixel-by pixel" 
comparisons of Ehi and Esfr. We first focus on stochasticity in 
the star formation and return to the comparison between E|^r and 
Ef^m Sec.EH 
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Figure 1. Overlays of the Ha, UV and HI images for UGC685. [A]Greyscales Ha (in 10 18 ergs s 1 per pixel of area 0.1225 arcsecond squared) , contours 
GALEX FUV image (from 0.0014 to 0.032 cps per pixel of area 2.25 arcsecond squared, in steps of 2). [B]Greyscales Ha, contours GMRT HI image (from 
17.5 to 1120 Jy/bmxm/s in steps of y/Z). [C]Greyscales GALEX FUV (in 10 -3 cps per pixel of area 2.25 arcsecond squared), contours GMRT HI image. 
Respective resolutions are, Ha: 1.9", FUV: 4", HI: 17" X 16". The length of the bold line in panel [A] is approximately 1 Kpc. 




Figure 2. [A]Esfr derived from Ha (Eg^^, empty squares) and FUV (Egpj^, filled circles) plotted against Shi, assumed to represent E gas , both 
axes being in log scale. The solid li ne represen t s the Kennicutt-Schmidt law with a slope of 1.4, and the dashed line represents the best fit Sch midt law fo r 
spiral galaxies only, both taken from Kennicutt ( 1998). The shaded region covers various estimates of the "threshold density" tabul ated in Kennicutil <1989h , 
[B]Disc-averaged values of aIK * ^SFR • ^ ne so ^ ' me ' s me 1 ; 1 une > an ^ me dashed line represents the relationship found bv lLee et el.l<20091) . [C]Ratio 

of global Egp^/Egp^as a function of Eg^^. The dashed line is the best fit straight line and has a slope of -0.63. The vertical dot-dashed line shows the 
approximate Esfr. value for our sample gala xies for which the SFR estimated assuming a Salpeter IMF will start deviating from the true SFR according to 
IPflamm-Altenburg. Weidner & Kroupal j2007l) . See the text for more details. 



3.1 Stochasticity in Star Formation 

R09 showed that from a comparison of the FUV and HI images, in 
star forming regions Efpp^ and Shi are related as 

log Egprf = (1-81 ± 0.05) log E H i - 4.70 ± 0.05 (3) 
By comparison with the canonical K-S law 

log Esfr = (1.4 ± 0.15) log E gas - 3.60 ± 0.14 (4) 

and noting that (i)Em is a strict lower limit to the total E gas and (ii) 
for a given FUV flux the inferred SFR decreases with decreasing 
metalicity, the robust conclusion that one can draw is that the star 
formation process in dwarf galaxies is significantly less efficient 
than that in big galaxies. R09 also showed (see their Fig. 6) that the 
data implied stochasticity and were best modelled as a stochastic 
power law with a variation of 50% in the coefficie nt (as opposed 
to the slope) of the power law. [Begum et al.| l j2006t ) had also high- 
lighted the stochasticity in the relation between Esfr and Ehi in 
dwarf galaxies. To properly characterize the star formation process, 
one would hence also need to know the average fraction of the gas 
that is participating in the star formation process. 

Figure[3]shows the fraction of pixels which are observed to be 
star forming (i.e. have a star formation rate of at least 3<j, where 
a is the rms in the UV image, in units of the star formation rate). 
The plot averages over 16 of the original sample of 23 galaxies, 7 



galaxies with relatively low GALEX exposure times have been ex- 
cluded. The dashed vertical lines indicate the rms level (after being 
translated from Esfr to Ehi using Eqn. |3j. For a given galaxy, 
if there was no stochasticity in the star formation, all points above 
the rms level (right of the corresponding dashed line) should have 
had observable star formation. As such, all points to the right of the 
rightmost dashed line can hence be regarded as giving a reliable 
fraction of gas that is participating in star formation. There are sev- 
eral points worth noting, viz. (1) all pixels with gas density greater 
than ~ IOMq/pc 2 participate in star formation. Interestingly, this 
number is identical to the thre shold density fo r star formation of 
~ 10 21 atoms/cm 2 proposed by |Skillmanl ( ll987l) , (2) the fraction of 
gas which participates in star formation decreases nearly linearly 
with decreasing Ehi (fsF = 0.96 logEm + 0.1). (3) even for a 
gas density logEm ~ —1.0, two orders of magnitude below the 
usually assumed threshold for star formation, at least 5% of the gas 
is observed to be forming stars. The average Ehi for the pixels with 
Ehi > 10M Q /pc 2 is ~ 17.3 M /pc 2 , and the average Esfr for 
these pixels is 3.5 x 10~ 3 M o /yr/kpc 2 . Thus even for the densest 
gas in dwarf galaxies, the star formation efficiency (i.e. EsFR/E gas ) 
is hence ~ 2.0~ 10 yr -1 , about a factor of tw o lower than the typical 
value for spiral galaxies dLerov et al .112008b . 

3.2 Massive star formation 

FUV emission is sensitive to the SFR of intermediate mass (M > 
3M0) relatively long lived (lifetime ~ 10 8 yr) stars. Ha emis- 
sion on the other hand traces the instantaneous SFR of massive 
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Figure 3. Plot of the fraction of gas that is "participating in star forma- 
tion" as a function of Shi- The plot averages over 16 galaxies, the dashed 
horizontal lines indicate the rms level of the UV images of the individual 
galaxies (after being translated from SsFR to Shi using Eqn. [3). All points 
to the right of the rightmost dashed line can be regarded as giving a reliable 
fraction of gas that is participating in star formation. The solid line is a fit 
to the "reliable" points. See text for more details. 
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Figure 4. Plots showing the binned 400pc resolution EsFR (hollow 
squares) and (half-filled circles) as a function of Shi- The Sgpjj 

data are for the 23 galaxies listed in R09, where as the data are for 

the 11 galaxies listed in Table [T] The average sensitivity levels for the two 
sets of data are indicated by the respective horizontal lines. The dashed and 
the dot-dashed lines show the Schmidt law fit to the Ha and FUV data 
respectively. The solid line is the canonical K-S law (Eqn.[4). The typical 
1 a scatter either above or below the mean in each bin is indicated by the 
vertical lines, in (1) the power law and (2) sensitivity limited parts for both 
tracers. Dashed lines are for Ha data, bold lines are for FUV data. Note that 
points with negative or zero FUV flux have been dropped from the plots. 



(M > 17M©) short lived (lifetime ~ 10 6 yr) stars. For our sample 
galaxies we show in Fig. [4] the Esfr as deduced from the FUV 
emission (Esfr'). as we H as Ha emission (Egjjy as a function 
of Ehi at a resolution of 400pc. (Note that pixels corresponding to 
gas not taking part in star formation, i.e. with negative or zero FUV 
flux, are not included in this plot) The best fit power laws to the Ha 
data is given by: 



logEfpR, = (1. 



± 0.04) log Ehi - 4.60 ± 0.05 (5) 



As can be seen the E SF r and E S fr points overlap within the 
scatter (indicated by the vertical line). Nonetheless as a comparison 
of Eqn. [5] and Eqn. [3] shows, there is a significant difference (~ 
2.7a, where a is the quadrature sum of the individual errors) in the 
slope of the two relationships, with the E|fr relation being steeper. 

Discrepancies between the SFR rates deduced be- 
tween these two tracers have been i nvestigated e arl ier b y 



several authors, including, for e.g. Meurer et aL d2009l); 



IPflamm-Altenburg. Weidner & Kroupal d2009h : iLee et el.l d2009l) 
A number of explanations for the two rates to diverge have been 



suggested, including (i) Stochastic paucity of high mass stars 
at low star formation rates. This would make Sj^j at low star 
formation rates lower than Esfr'' . For example ILee et el ] d2009l) 
show that for star formation rates lower than ~ lO _2 M0/yr, 
the Ha emission systematically under predicts the true SFR. 
(ii) Non uniform star formation rates. For example if the star 
formation is bursty, then a few million years after the burst all 
the OB stars would have died and the Ha emission would once 
again systematically underestimate the true SFR. (iii) leakage of 
ionizing photons, either out of the galaxy, or into a more diffuse 
region of the ISM, where the resultin g Ha emission has to o low a 
surface brightness to be detected (e.g. Mel ena et alj d2009h v Once 
again, this would result in the Ha emission underes timating the 
true S FR. (iv) variations in the IMF. For example, iMeurer et al.1 
(2009) identify correlations between Eg^/Eg^with global 
galaxy parameters like the luminosity, rotational velocity and 
dynamical mass, and argue that this implies an IMF that varies 
with environment. Weidner & Kroupa d2005h present a model in 
which the underlying IMF is universal, but a dependence of the 
most massive star formed in a cluster on the mass of the cluster 
leads to the total stellar population having a steeper IMF than the 
canonical one. (v) Dust extinction. Since dust extinction is more 
at the shorter wavelengths, under correction for dust would lead 
to the FUV emission underestimating the true SFR. Note that in 
most of the above scenarios the Ha emission would under predict 
the true star formation rate. One would expect the E|fr to exceed 
Eg^^ (as observed for about half of our sample) only if (i) the 
IMF is more top heavy than assumed, or (ii) the dust extinction has 
been underestimated. 

To explore this issue further, we show in Figure [5] pixel by 
pixel correlations of E|^ R /E| F ^ / (both at 150 pc resolution) with 
Esfr and Em(at 400 pc resolution). In each panel the hollow 
squares are for those galaxies for which the global E|^r is greater 
than the global EfpR , while the filled circles are for those galax- 
ies for which the global Eg^is less than the global Efpji , From 
the first panel, one can clearly see that the anti correlation between 
EgiTR/EgpREUid Esfr seen on global scales continues even on 
scales as small as 150pc. The right axis of the panels is the amount 
of differential dust obscuration required to bring the two SFR es- 
timators into agreement. From Figure [5]one can see that bringing 
the two SFR estimators into agreement at the lowest star forma- 
tion rates requires the dust obscuration to be more at Ha than at 
FUV, which is physically implausible. It is more likely that one of 
the several mechanisms discussed above for suppressing the Ha 
flux at low star formation rates is operative. At high star forma- 
tion rates, where E|^r > Esfr', the average Nhi/Av required 
to bring the two estimators into agreement is 8, i.e. the gas should 
be ab out twice as dust rich as the SMC (for which Nhi/Av is 16.3 
from lBouchet et alj dl985h ). If one assumes that these regions have 
substantial molecular gas, and that the ga l axies follow the L-Z re- 
lation for dwarfs (e.g. Ekta & Chengalur (2010)), and that dust is 
proportional to metalicity, then the required molecular gas densi- 
ties to bring the gas to dust ratio to the same value as the SMC is 
Eh 2 > lO 2 M0/pc -2 , similar to the peak densities in the center 
of spirals, which again seems unlikely. In summary it does not ap- 
pear that dust extinction is the primary cause of the disagreement 
between Eg^R and E| F r / at the high Esfr end. 

In terms of direct observables, panel [B] shows that for 
the same amount of FUV emission, galaxies with lower global 
Egpji^SFR are under producing Ha emission. This could either 
be because the galaxies have a fading st arburst or becau se the 
galaxies are not producing high mass stars. ILee et el.l J2009II3) find 
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Figure 5. The ratio of Eg^ and Sf^f (at scales of 150 pc) plotted as a function of Eg^L, (Panel[A]), E|^ (Panel[B]) and E HJ (Panel[C]). The S H i is 
computed from 400pc resolution images. The solid dots are for for galaxies for which the global EjjjjS' > ^gpj^i an d the hollow square are for the galaxies 
with global Sf p a R > Efjjig' . Means and the en'ors on the mean for each bin are shown. The axes on the right show the amount of differential dust extinction 
required to bring the two SFR estimators into agreement. When binning in Eg^ all pixels with E|^ more than 3 times the rms value (of the corresponding 
galaxy) are considered (regardless of the value of Eg^^), similarly, when binning in Eg^ 7 (and binning in Ehi) all pixels with EgS^ more than 3 times 
the rms value are considered (regardless of the value of E|^y. In Panel[A] the dotted rectangular region marks the area covered in the similar global plot 
FigureffJC]. The dashed line in it is the best fit straight line from the global plot FigureEtC]. 



that the frequency and amplitude of star bursts in dwarfs make 
the former explanation unlikely. However, a more detailed calcula- 
tion, and observations of a larger sample would be needed to prop- 
erly settle this issue. The most striking feature of the plots how- 
ever is in panel [C], which shows that galaxies with lower global 
Egpja/EgFR do not have gas with column density > lOM0/yr. 
The most straight forward interpretation of this is that massive star 
formation is more likely to happen in gas with high column den- 
sities. Inde ed, star formation model s have supported such a corre- 
lation (e.g. iKrumholz etaD fcOldV ). While the linear scales that 
the models refer to are much smaller than those that we are deal- 
ing with here, such a correlation is likely given that high density 
star forming regions are more likely to occur in regions where the 
overall gas density is higher. 



4 CONCLUSIONS 

We find a clear stochasticity between the Ehi and Esfr. All gas 
with Ehi > lOM© pc -2 has associated star formation. While 
the fraction of star forming gas decreases with decreasing Ehi 
there is no sharp "threshold" below which star formation is com- 
pletely quenched. We also find that galaxies for which globally 
^SPR < ^sfr are marked by not having high HI column den- 
sity (i.e. Ehi > lOM© pc~ 2 ) gas. This is consistent with models 
in which formation of high mass stars preferentially happens in re- 
gions with high gas column density. 
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